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Abstract 

We show that the effective field theory of low energy modes in dense QCD has 
positive Euclidean path integral measure. The complexity of the measure of 
QCD at finite chemical potential can be ascribed to modes which are irrelevant 
to the dynamics at sufficiently high density. Rigorous inequalities follow at 
asymptotic density. Lattice simulation of dense QCD should be possible using 
the quark determinant calculated in the effective theory. 
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Quantum chromodynamics (QCD) with a non-zero chemical potential has a complex 
measure, which has thus far precluded lattice simulations [Q. Recent analytical work in 
color superconductivity has demonstrated a rich phase structure at high density, and 
stimulated interest in QCD at non-zero baryon density Several experiments have been pro- 
posed to probe matter at density of a few times nuclear matter density [[| . Even rudimentary 
information about the behavior of dense matter would be useful to the experimental pro- 
gram, as well as to the study of compact astrophysical objects such as neutron stars. In 
this paper, we will show that QCD near a Fermi surface has positive, semi-definite measure. 
The contribution of the remaining modes far from the Fermi surface can be systematically 
expanded, using a high density effective theory previously introduced by one of us Q. This 
effective theory is sufficient to study phenomena like color superconductivity, although quan- 
tities like the equation of state are presumably largely determined by dynamics deep in the 
Fermi sea. 

Let us recall why the measure of dense QCD is complex in Euclidean space. We use the 
following analytic continuation of the Dirac lagrangian to Euclidean space: 

x -> -ix%, Xi -> x l E ; 70 -> 7I, 7* -> vf E . (1) 

The Euclidean gamma matrices satisfy 

7^ = 7£ , {l^Y E } = 2^. (2) 

The Dirac-conjugated field, if) = ip^j , is mapped into a field, still denoted as if), which is 
independent of ip and transforms as if)' under SO (A). Then, the grand canonical partition 
function for QCD is 

Z(p) = J dA^ det (M) e" 5 ^, (3) 

where S(A fi ) is the positive semi-definite gauge action, and the Dirac operator 

M = 1e D e + /i 7 |, (4) 

where D E = Oe+^Ae is the analytic continuation of the covariant derivative. The Hermitian 
conjugate of the Dirac operator is 

M^ = - lE D E + mE . (5) 

The first term in (|J) is anti-Hermitian, while the second is Hermitian, hence the generally 
complex eigenvalues. When \i = 0, the eigenvalues are purely imaginary, but come in 
conjugate pairs (A, A*) [§], so the resulting determinant is real and positive semi-definite: 

detM = I]A*A>0 . (6) 

In what follows we investigate the positivity properties of an effective theory describing 
only modes near the Fermi surface. A system of degenerate quarks with a net baryon number 
asymmetry is described by the QCD Lagrangian density with a chemical potential fi, 

£ QCD = - \F a , u F a ^ + ^ 7o ^, (7) 
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where the covariant derivative = <9 M + iA^ and we neglect the quark mass for simplicity. 
The energy spectrum of (free) quarks is given by an eigenvalue equation, 

(5-p + /i)^± = E±ij;±, (8) 

where a = 707 and ip± denote the energy eigenfunctions with eigenvalues E± = —\i ± |p|, 
respectively. At low energy E < fi, the states ■?/>+ near the Fermi surface, \p\ ~ //, are easily 
excited but ■?/>-, which correspond to the states in the Dirac SGcL, £1X6 completely decoupled 
due to the presence of the energy gap \x provided by the Fermi sea. Therefore the appropriate 
degrees of freedom at low energy consist of gluons and ip + only. 

Now, we wish to construct an effective theory describing the dynamics of ip + by in- 
tegrating out modes whose energy is greater than /1. Consider a quark near the Fermi 
surface, whose momentum is close to fj,v F . Without loss of generality, we may decompose 
the momentum of a quark into a Fermi momentum and a residual momentum as 

Vp. = /^V + V' ( 9 ) 
where v M = (0, v F ). Since the quark energy is given as 



E = -n + y/m + + 11, (10) 

the residual momentum should satisfy (In + fi) 2 + l\< 4/i 2 with In = v F l ■ v F and l± — I — In. 

To describe the small excitations of the quark with Fermi momentum, /j,v f , we decompose 
the quark fields as 

^( x ) = e i ^ s [i/; + (vp,x)+^(v F ,x)} , (11) 

where 

il>±(vF,x) = P±(v F )e- i '**--*il>(x) with P ± (v F ) = l±a 2 VF . (12) 

The quark Lagrangian in Eq. (|7|) then becomes 

ip (ip+ /i7°) = i>+(v F , x)i^ D^i/j^vp, x) + ^_({7 F ,x)7° \ 2fi + iD^j ip_(v F ,x) 

+ tjj-(v F ,x)ip ± i/j+(v F ,x) + h.c. (13) 

where 7J = (7 , v F v F ■ 7), 7^ = 7 M -7f, D\\ = V^D^ with W = (l,v F ), = (l,-v F ), and 

IP. -"/v 

At low energy, we integrate out all the "fast" modes ip- and derive the low energy effective 
Lagrangian by matching all the one-light particle irreducible amplitudes containing gluons 
and in loop expansion. The effects of fast modes will appear in the quantum corrections 
to the couplings of low energy interactions. At tree-level, the matching is equivalent to 
eliminating ip_ in terms of equations of motion: 

• • 00 / ■ j-y \ n 

x) = -whsf^, *) = E (-VJ «'+ ( *- *>■ (14) 



3 



Therefore, the tree-level Lagrangian for ip + becomes 

C° eS = - ^+7°(#l) V+ + • • • , (15) 

where the ellipsis denotes terms with higher derivatives. 

Consider the first term in our effective Lagrangian, which when continued to Euclidean 
space yields the operator 

M eft = 7 f ■ D{A) . (16) 

M e f t is anti-Hermitian and it anti-commutes with 75, so it leads to a positive semi-definite 
determinant. However, note that the Dirac operator is not well defined in the space of 
i/} + (vf,x) (for fixed vf), since it maps ip + (vF,x) into if)+(— vp, x): 

ilp\\P + ^ = P-%Jp\\i\) . (17) 

Since P-{vf) = P+(~vf), ip */>+(vf, x ) are i J +(—'VF, x) modes, or fluctuations of a quark 
with momentum —fivp. 

We can demonstrate the necessity of including both ip + (vF,x) and ip + (—VF,x) modes in 
our effective theory by considering charge conservation in a world with only +vp quarks. 
The divergence of the quark current at one loop is 

(d, J a »(v F , x)) =g s f ^e-^Iif u {p)A b ({-p) , (18) 

— * , 

where A\\ = (A , VpVp • A) and 11^, is the vacuum polarization tensor in the effective theory 
given as || 

TV ,u ( x • 2 A P-VPV»V U nQ , 

Kb (P) = W $ab „ : . „ z; ■ 19 

p ■ V + lep ■ v p 

The polarization tensor has to be transverse to maintain gauge invariance. We find that 
if we have both fields t/) + (vp,x) and vf, x) the current is conserved and the gauge 
symmetry is not anomalous: 

(d^(v F ,x) + d^(-Vp,x)) = 0. (20) 

Therefore, we need to introduce quark fields with opposite momenta. The Dirac operator is 
well defined on this larger space. 

This anomaly can be understood in terms of spectral flow, since the Fermi surface is 
(in a certain sense) not gauge- invariant. Under a gauge transformation, U(x) = e iq ' x , the 
Hamiltonian changes and the energy spectrum of free modes of residual momentum I shifts 
to E = I ■ vp + q ■ vp. Quarks near the Fermi surface with vp ■ q > flow out of the Fermi 
sea, creating charge. This charge creation is compensated by quarks with opposite Vp; their 
energy decreases and they flow into the Fermi sea. However, unless modes with opposite 
velocities (i.e. both sides of the Fermi sphere) are included, charge is not conserved. 

Thus far we have considered the quark velocity as a parameter labelling different sectors 
of the quark field. This is similar to the approach of heavy quark effective theory (HQET) || , 
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in which the velocity of the heavy charm or bottom quark is almost conserved due to 
the hierarchy of scales between the heavy quark mass and the QCD scale. However, this 
approach contains an ambiguity often referred to as "reparameterization invariance" , related 
to the non-uniqueness of the decomposition (Qj of quark momenta into a large and residual 
component. In the dense QCD case, two ip(vp, x) modes whose values of vf are not very 
different may actually represent the same degrees of freedom of the original quark field. In 
what follows we give a different formulation which describes all velocity modes of the quark 
field, and is suitable for defining the quasiparticle determinant. 

First, a more precise definition of the breakup of the quark field into Fermi surface modes. 
Using the momentum operator in a position eigenstate basis: p = —id, we construct the 
Fermi velocity operator: 

—id 
V ~ V /Zy2 

which is Hermitian, and a unit vector. 

Using the velocity operator, we define the projection operators P± as before and break 
up the quark field as, if>(x) = if)+(x) + if)-(x), with if)± = P±if). By leaving v as an operator 
we can work in coordinate space without introducing the HQET-inspired velocity Fourier 
transform which introduces vp as a parameter. If we expand the quark field in the eigenstates 
of the velocity operators, we recover the previous formalism with all Fermi velocities summed 
up. 

The leading low-energy part of the quark action is given by 

As before, we define the fields if> + to absorb the large Fermi momentum: 

ifj + (x) = e- iflM P + (v)iP{x). (23) 

Let us denote the eigenvalue v obtained by acting on the field if) (which has momentum 
of order p) as v\ (or v "large"), whereas eigenvalues obtained by acting on the effective 
field theory modes ip + are denoted v r (or v "residual"). If the original quark mode had 
momentum p with \p\ > p (i.e. was a particle), then vi and v r are parallel, whereas if \p\ < p 
(as for a hole) then v r and vi are anti-parallel. In the first case, we have P+{vi) = P + (v r ) 
whereas in the second case P+(vi) = P-{v r ). Thus, the residual modes can satisfy either 
of P±(v r )ip+ = depending on whether the original ip mode from which it was derived 
was a particle or a hole. In fact, ip + modes can also satisfy either of P±{vi)ip+ = V ; + since 
they can originate from if) modes with momentum ~ +pv as well as —pv (both are present 
in the original measure: Dip Dip). So, the functional measure for ip + modes contains all 
possible spinor functions - the only restriction is on the momenta: |Zo|, |^| < A, where A is 
the cutoff. 

In light of the ambiguity between vi and v r , the equation if> = e +l ^ x ' v if> + must be modified 

to 

if) = exp (+ipx ■ v a ■ v)ip + = exp (+ipx • v r a • v r ) ip + , (24) 
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(21) 



(22) 



where the factor of a • v r corrects the sign in the momentum shift if v r and v\ are anti- 
parallel. In general, any expression with two powers of v is unaffected by this ambiguity. 
For notational simplicity we define a local operator 

oc l x^ d 2 

X = ^x-va-v = HyTQjQj. (25) 

Taking this into account, we obtain the following action: 

jr, + = $ + e- iX (i$-4 + Wo)e+ iX ij; + . (26) 

We treat the A" term separately from + //70 since the former does not commute with X, 
while the latter does. Continuing to Euclidean space, and using the identity P_ 7aj P + = 7IP+, 
we obtain 

C + = V5 +7 f (fl" + iA$) ip + , (27) 

where 

A% = e~ iX A» e +lX , (28) 

and all 7 matrices are Euclidean. The term containing A cannot be fully simplified because 
[v, A] 7^ 0. Physically, this is because the gauge field carries momentum and can deflect the 
quark velocity. The redefined tf) + modes are functions only of the residual momenta 1, and 
the exponential factors in the A term reflect the fact that the gluon originally couples to the 
quark field ip, not the residual mode ip+- 

The kinetic term in ( p7j ) can be simplified to 

= rf& (29) 

since v-d v -7 = d-'y . The action (^) is the most general dimension 4 term with the rotational, 
gauge invariance lil and projection properties appropriate to quark quasiparticles. Therefore, 
it is a general consequence of any Fermi liquid description of quark-like excitations. 

The operator in (^7|) is anti-Hermitian and leads to a positive, semi-definite determinant 
since it anti-commutes with 75. The corrections given in (|1^) are all Hermitian, so higher 
orders in the 1/fi expansion may re-introduce complexity. The structure of the leading 
term plus corrections is anti-Hermitian plus Hermitian, just as in the original QCD Dirac 
Lagrangian with chemical potential. The leading terms in the effective action for gluons 
(these terms are generated when we match our effective theory, with energy cutoff A, to 
QCD) also contribute only real, positive terms to the partition function: 

r ( 1 M 2 \ 

S eS (A) = J d 4 x E l^F^ + — E A%A% \ > 0, (30) 



1 If we simultaneously gauge transform A + and V+ in (P^l) the result is invariant. There is a 
simple relation between the gauge transform of the + fields and that of the original fields: U+(x) = 
U(x)e lX . Of course, the momentum-space support of the + gauge transform must be limited to 
modes less than the cutoff A. 
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where A± — A — A\\ and the Debye screening mass is M = ^jNf/(2ir 2 )g s fi . 

Matching of hard gluon effects also leads to four-quark operators in the effective theory. 
At asymptotic density, we can neglect these operators, since forward scattering dominates 
Cooper pairing interactions (due to Landau damping 0). However, at lower densities hard 
operators may be important. Matching effects due to hard gluon exchange still lead to a 
positive action for attractive channels, since they arise from quasiparticle-gluon interactions 
which are originally positive!. Only interactions involving virtual anti-quarks lead to non- 
positive interactions, and these are always suppressed by powers of \i. 

Positivity of the measure allows for rigorous QCD inequalities at asymptotic density. For 
example, inequalities among masses of bound states can be obtained using bounds on bare 
quasiparticle propagators. One subtlety that arises is that a quark mass term does not lead to 
a quasiparticle gap (the mass term just shifts the Fermi surface). Hence, for technical reasons 
the proof of non-breaking of vector symmetries [7| must be modified. (Naive application of 
the Vafa-Witten theorem would preclude the breaking of baryon number that is observed 
in the color-flavor-locked (CFL) phase ||.) A quasiparticle gap can be inserted by hand 
to regulate the bare propagator, but it will explicitly violate baryon number. However, 
following the logic of the Vafa-Witten proof, any symmetries which are preserved by the 
regulator gap cannot be broken spontaneously. One can, for example, still conclude that 
isospin symmetry is never spontaneously broken. In the case of three flavors, one can use 
the CFL gap as a regulator to show rigorously that none of the symmetries of the CFL 
phase are broken at asymptotic density. On the other hand, by applying anomaly matching 
conditions 0, we can prove that the axial symmetries are broken. We therefore conclude 
that the CFL phase is the true ground state for three light flavors at asymptotic density. 

It may be possible to simulate dense QCD using the effective field theory determinant in 
place of the usual quark determinant. We know that this is a good approximation at very 
high density, and it should remain a good approximation as long as Aqcd, the characteris- 
tic scale of the dominant physics, is smaller than //. At lower densities the weak coupling 
approximation no longer holds, so (13) is no longer a good guide to the higher order correc- 



tions, which are only constrained by symmetry and projection requirements. However, we 
can use naive dimensional analysis [IU| to estimate the coefficients of the higher dimension 



operators. 

An estimate of the size of corrections to the determinant from higher orders can be 
obtained by considering the Euclidean relation 

detM = e™ nM = e" eoV , (31) 

where €q is the vacuum energy density and V the volume of the system. Corrections to 
the vacuum energy can be estimated using conventional diagrammatic methods and naive 
dimensional analysis. We find that corrections are roughly suppressed by powers of 



2 A simple way to study the positivity of four-quark operators is to replace them by a vector field 
with trivial quadratic term V 2 which couples to quarks like the original gluon: V^-tp^/^ip. Completing 
the square, we see that the resulting path integral is positive if the four-quark interactions are 
attractive. 
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To realize the effective theory ( p7|) directly on the lattice, one can replace the plaquette 
U ntn+fl ~ 1 + iaA^in) by U + = e~ tX Ue +tX in the fermion action, but not in the gauge 
action. In effect, one computes the fermion determinant in the usual way, but as a function 
of U + . The momenta of the quasiparticle ip + modes is simply the residual momenta 1, which 
is unrestricted except that its magnitude must be small compared to \i. One can impose 
this condition on ip + by choosing a lattice spacing ad e t >> A^ 1 in the quark determinant. 
However, a challenging aspect of this determinant is that it is computed on the geometry of 
a spherical shell rather than a ball: 

d 3 p = dp p 2 dQ = dl(n + I) 2 dQ dH = dl I 2 dVt . (32) 

To obtain the e ±lX operators, one must first realize (and presumably diagonalize) the 



velocity operator (£L|). Since the momentum operator has a simple representation in coor- 
dinate space, the most challenging aspect of v is the normalization factor (the square root 
of the Laplacian in (|21f)), which is non-local. (A similar problem arises in lattice models of 
chiral fermions.) One should probably investigate a simpler method for directly enforcing 
the normalization condition on v. Note, however, that v and X are independent of the gauge 
field, so they need be computed only once for each lattice. 

Probably the most direct way to utilize our results on the lattice is as follows. Imagine 
coupling dense quark matter to a background gauge field A whose magnitude and derivatives 
are characterized by a scale A << \x. The leading part of the low-energy effective theory 
(describing only Fermi surface modes) has a real and positive determinant. The quark 
determinant (or equivalently, its logarithm which is the effective action) can be expanded 
in powers of with the leading term real and positive. Since the determinant is a 
functional of the gauge field and its derivatives, the expansion will effectively be in powers 
of A over fi. This means that the ordinary lattice determinant det^-D^ + ^74] computed in 
such backgrounds should be real and positive to leading order in Physically, the low- 
momentum gauge fields cannot excite modes deep within the Fermi sea (such as anti-quarks) 
which lead to complex contributions to the determinant. If we further take A >> Aqcd, 
then the neglected effects of higher momentum gluon modes are suppressed by powers of 
a, (A) « I. 

In order to restrict ourselves to gauge fields satisfying A << fi we need to couple the 
quark determinant living on a lattice with spacing a ~ to a set of gauge fields living 
on a coarser lattice of spacing a' ~ 1/A. Ordinary calculations of the dense matter quark 
determinant on lattices with a = a' will not yield a real, positive result since the sub-leading 
terms in the expansion of the determinant are as large as the leading term. 
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